Noise and Inertia-Induced Inhomogeneity in the Distribution of Small Particles in Fluid Flows 
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The dynamics of small spherical neutrally buoyant particulate impurities immersed in a two-dimensional fluid 
flow are known to lead to particle accumulation in the regions of the flow in which rotation dominates over shear, 
provided that the Stokes number of the particles is sufficiently small. If the flow is viewed as a Hamiltonian 
dynamical system, it can be seen that the accumulations occur in the nonchaotic parts of the phase space: the 
Kolmogorov-Arnold-Moser tori. This has suggested a generalization of these dynamics to Hamiltonian maps, 
dubbed a bailout embedding. In this paper we use a bailout embedding of the standard map to mimic the 
dynamics of impurities subject not only to drag but also to fluctuating forces modelled as white noise. We 
find that the generation of inhomogeneities associated with the separation of particle from fluid trajectories is 
enhanced by the presence of noise, so that they appear in much broader ranges of the Stokes number than those 
allowing spontaneous separation. 

PACS numbers: 47.52.+j, 05.45.Gg, 45.20.Jj 



I. INTRODUCTION 

Impurities suspended in a fluid flow are frequently observed 
to be distributed inhomogeneously. Even in very chaotic 
flows, particulate impurities arrange themselves in extraordi- 
narily structured distributions, in apparent contradiction to the 
high mixing efficiency expected from the characteristics of the 
basic flow. To give just one example, in the particular instance 
of geophysical fluids, the filamentary structure, or patchiness, 
often displayed by plankton populations in the oceans is apuz- 
zling problem currently under intense investigation [|l], ^ ||]. 
Several mechanisms to produce this type of inhomogeneity 
have been studied, and include both dynamical aspects of the 
flow as well as the reactive properties of the considered im- 
purities. The basic idea in these mechanisms is that the par- 
ticle loss due either to the flow — in open flows — or to the 
chemical or population dynamics of the particles — in closed 
flows — is minimized on some manifolds associated with the 
hyperbolic character of the flow H ^jl. In this paper we 
explore an alternative purely dynamical mechanism for inho- 
mogeneity with non-reactive particles in bounded flows. We 
show that particle inertial effects combined with fluctuating 
forces are capable of producing inhomogeneity even in cases 
in which the impurity and fluid densities match exactly. 

When impurities have a different density to the fluid, it is 
intuitively clear that they will be expelled from rapidly rotat- 
ing regions of the flow — for heavy particles — or attracted 
to the center of those regions — for light particles — because 
of centrifugal effects [0]. However, it was recently demon- 
strated that neutrally buoyant particles also tend to settle in the 
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rotation-dominated regions of a flow, but because of a more 
subtle mechanism involving the separation between the fluid 
and particle trajectories that can occur in the opposite regions, 
i.e., in the areas of the flow dominated by shear [jsj. How- 
ever, this mechanism is only relevant when the particle Stokes 
number is smaller than the eigenvalues of the Jacobian ma- 
trix of the flow, a condition that may not be fulfilled in some 
physically interesting situations. We show here by means of 
a minimal model that the latter condition may be relaxed if 
a small amount of noise be added to the forces acting on the 
impurity. 

Our approach is qualitative, in the sense that instead of 
considering a specific flow and the precise particle dynamics 
induced by it, we describe the system with an iterative map 
whose evolution contains the basic features of both the fluid 
flow and the particle dynamics: flow volume preservation, to- 
gether with particle separation at the hyperbolic regions. The 
reason for moving to a discrete system is that in the map the 
phenomena that we describe may be understood more intu- 
itively, while translating the results back to the flow case is 
immediate. The strategy of understanding fluid-dynamical 
phenomena by using iterated maps amenable to the powerful 
artillery of dynamical-systems theory has been successfully 
applied on several different occasions. For example, the struc- 
tures of the chaotic advection induced by time-periodic three- 
dimensional incompressible flows were predicted by studying 
the qualitatively equivalent dynamics of three-dimensional 
volume-preserving maps ^|. Later, these structures were 
confirmed in realistic flows [jn], 12 13|]. Other examples are 
the treatment of the propagation of combustion fronts in lami- 
nar flows by a qualitative map approach [|l4|, and the descrip- 
tion of the formation of plankton population structures due to 
inhomogeneities of the nutrient sources [[T^. Remarkably, in 
the instance we discuss in this paper, the procedure is also use- 
ful from the point of view of dynamical-systems theory, as it 
has suggested a new technique — bailout embedding — for 
the control of Hamiltonian chaos 



The plan of the paper is as follows. First, we briefly review 
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the classical model for the forces acting on a small spheri- 
cal tracer moving relative to the fluid in which it is immersed 
(Section II). Concentrating on the case of a neutrally buoyant 
impurity, we trace the construction of a minimal model that 
makes evident the separation of particle and fluid trajectories 
in the regions in which the flow presents strong shear (Sec- 
tion |l|). On the basis of this model, we make a generalization 
that allows us to build a discrete mapping that represents the 
Lagrangian evolution of the fluid parcels as well as the dynam- 
ics of the particle (Section While this map also displays 
particle-fluid separation when the parameter 7, equivalent to 
the Stokes number, is relatively small, we show in Section ^ 
that a small amount of noise, added to the dynamics of the 
particle to separate it continually from the flow, enhances the 
impact of the hyperbolic regions far beyond the values of 7 
required for separation. Conclusions are to be found in Sec- 
tion Ivil 



II. MAXEY-RILEY EQUATIONS 

The equation of motion for a small, spherical tracer in an 
incompressible fluid we term the Maxey-Riley equation 
O, Eq], which may be written as 
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Here v represents the velocity of the particle, u that of the 
fluid, pp the density of the particle, pf, the density of the fluid 
it displaces, v, the kinematic viscosity of the fluid, a, the ra- 
dius of the particle, and g, gravity. The derivative Du / Dt is 
along the path of a fluid element 



Du du , 

— = - + («. V)n, 



(2) 



whereas the derivative du/dt is taken along the trajectory of 
the particle 



(3) 



du du , , 

- = - + V)n. 

The terms on the right of Eq. (|l|) represent respectively the 
force exerted by the undisturbed flow on the particle, buoy- 
ancy. Stokes drag, the added mass due to the boundary layer 
of fluid moving with the particle jl^ pO| ], and the Basset- 
Boussinesq force [ pTj , ^ | that depends on the history of 
the relative accelerations of particle and fluid. The terms in 
a^V^tt are the Faxen corrections. The Maxey-Riley 
equation is derived under the assumptions that the particle ra- 
dius and its Reynolds number are small, as are the velocity 
gradients around the particle. 



III. MINIMAL MODEL 

First let us consider a minimal model for a neutrally buoy- 
ant particle. For this we set Pp = pj in Eq. (|l]). We consider 
the Faxen corrections and the Basset-Boussinesq term to be 
negligible. We now rescale space, time, and velocity by scale 
factors L,T = L/U, and U, to arrive at the expression 

dv Du „ 1 , , 1 f dv Du\ 

dT-i^-^^ (^-«)-2U"i5rj' 

where St is the particle Stokes number St = 2a^U / {%vL) = 
2/9 (a/L)^Re/, Re/ being the fluid Reynolds number The 
assumptions involved in deriving Eq. (|l]) require that Si <C 1 
in Eq. (|). 

If we substitute the expressions for the derivatives in Eqs. 
(||) and (^) into Eq. (^, we obtain 

^{v-u)^-{(v-u)-V)u-l St-' (v-u). (5) 
dt 3 



We may then write A = v — u, whence 
dA / 2 1 \ 



(6) 



where J is the Jacobian matrix — we now concentrate on two- 
dimensional flows U = {Ux, Uy) 



J 



If we diagonalize the matrix we obtain 

dAp _ ( A-2/3St"^ 

-A-2/3St-i 



dt 







(7) 



Ad, (8) 



so if Re{X) > 2/3 St^^, Ad may grow exponentially. Now 
A satisfies det( J — AI) — 0, so A-^ — trJ + detJ ~ 0. Since 
the flow is incompressible, dxUx + dyUy — tr J — 0, thence 
— A^ = dctJ. Given squared vorticity uj'^ = {dxUy — dyUxY, 
and squared strain = + s^, where the normal component 
is si = dxUx — dyUy and the shear component is S2 = dyUx + 
dxUy, we may write 

Q = A^ = -detJ = {s^ - u?)l\, (9) 

where Q is the Okubo-Weiss parameter @, |5p. If Q > 0, 
A^ > 0, and A is real, so deformation dominates, as around 
hyperbolic points, whereas if Q < 0, A^ < 0, and A is com- 
plex, so rotation dominates, as near elliptic points. Equation 
(jg) together with dx /dt = A + u defines a dissipative dy- 
namical system 



d^/dt = F(^) 



(10) 



with constant divergence V .F = —4/3 St^ in the four di- 
mensional phase space ^ = (a;,y,ylj^,ylj,), so that while small 
values of St allow for large values of the divergence, large val- 
ues of St force the divergence to be small. The Stokes number 
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is the relaxation time of the particle back onto the fluid trajec- 
tories compared to the time scale of the flow — with larger St, 
the particle has more independence from the fluid flow. From 
Eq. (||), about areas of the flow near to hyperbolic stagnation 
points with Q > 4/9 St^^, particle and flow trajectories sep- 
arate exponentially. 
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FIG. 1 : (a) Contour plot illustrating magnitude of Q — lighter is 
higher Q — for the time-independent model Eq. ( [ll| ) (the flow is on 
a torus), (b) The separation of a neutrally buoyant particle trajectory 
(red line) from the flow (blue line) in regions of high Q allows the 
particle to wander between cells, (c) After a complicated excursion, 
a particle (red line) eventually settles in a zone of low Q of the flow 
(blue line), (d) The velocity difference — Ux between the particle 
and the flow against time. 

To illustrate the effects of S< and Q on the dynamics of a 
neutrally buoyant particle, let us consider the simple incom- 
pressible two-dimensional model flow defined by the stream 
function 



tp{x,y,t) — ^cos(x + B sin cut) cosy. 



(11) 



The equations of motion for an element of the fluid will then 
be i = dyip, y = —dxip- "0 has the role of a Hamiltonian 
for the dynamics of such an element, with x and y playing the 
parts of the conjugate coordinate and momentum pair. Let us 
first consider the simplest case, for which we suppress time 
dependence by setting B = 0. Thence ip should be a constant 
of motion, which implies that real fluid elements follow tra- 
jectories that are level curves of ijj. In Fig. |l](a) are depicted 
contours of Q. Notice that the high values of Q are around the 
hyperbolic points, while negative Q coincides with the cen- 
tres of vortices — elliptic points — in the flow. Figure [l| (top 
right) shows the trajectory of a neutrally buoyant particle start- 
ing from a point on a fluid trajectory within the central vortex. 



but with a small velocity mismatch with the flow. This mis- 
match is amplified in the vicinity of the hyperbolic stagnation 
points for which Q is larger than 4/9 Si~^ to the extent that 
the particle leaves the central vortex for one of its neighbors, 
a trip that is not allowed to a fluid parcel. In the end a par- 
ticle settles on a trajectory that does not visit regions of high 
Q, proper for a fluid parcel. While this effect is already seen 
in Fig. |I|(b), it is more dramatically pictured in the trajectory 
shown in Fig. 0(c), in which the particle performs a long and 
complicated excursion wandering between different vortices 
before it settles in a region of low Q of one of them. To il- 
lustrate the divergence of particle and fluid trajectories, and 
the fact that particle and fluid finally arrive at an accord, in 
Fig. |l](d) we display the difference between the particle veloc- 
ity and the fluid velocity at the site of the particle against time 
for this case. Notice that this difference seems negligible at 
time zero, and that it also convergences to zero at long times, 
but during the interval in which the excursion takes place it 
fluctuates wildly. 

Even more interesting is the case of time-dependent flow: 
_B 7^ in our model. As in a typical Hamiltonian system, as- 
sociated with the original hyperbolic stagnation points, there 
are regions of the phase space — here real space — domi- 
nated by chaotic trajectories. Trajectories of this kind, stro- 
boscopically sampled at the frequency of the flow, are re- 
produced in Fig. g. Such trajectories visit a large region of 
the space, which includes the original hyperbolic stagnation 
points and their vicinities where Q is large. Excluded from 
the reach of such a chaotic trajectory remain areas where 
the dynamics is regular; the so-called KAM (Kolmogorov- 
Arnold-Moser) tori. In our model these lie in the regions 
where Q < 4/9St~^. Now a neutrally buoyant particle try- 
ing to follow a chaotic flow pathline would eventually reach 
the highly hyperbolic regions of the flow. This makes likely 
its separation and departure from such a pathline, in search 
of another pathline to which to converge. However, conver- 
gence will only be achieved if the pathline never crosses areas 
of high Q. Figure ^ demonstrates this phenomenon: particles 
were released in the chaotic zone with a small velocity mis- 
match. Such a particle follows the flow, until, coming upon 
a region of sufficiently high Q, it is thrown out of that flow 
pathline onto a long excursion that finally ends up in a regular 
region of the flow on a KAM torus. The regular regions of 
the flow then constitute attractors of the dissipative dynami- 
cal system Eq. ( p^ ) that describes the behaviour of a neutrally 
buoyant particle. 



IV. DISCRETE DYNAMICS DESCRIPTION 

An examination of the dynamical system defining our min- 
imal model for the behavior of neutrally buoyant particles 
shows that it is composed of some dynamics within some 
other larger set of dynamics. Eq. (^ can be seen as an equa- 
tion for a variable that would define another equation of mo- 
tion, when the solution of the former be zero. Recall that the 
inner dynamics is that representing the Lagrangian trajecto- 
ries of the flow, and that in two dimensions, it is a Hamil- 
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FIG. 2: Poincare sections of trajectories in the time-dependent flow 
of Eq. ([h]). From top left to bottom riglit are siiown (red dots) four 
increasingly chaotic examples of the flow, and (magenta crosses), 
the trajectories of neutrally buoyant particles in the flows that in each 
case finally end up on a KAM torus within the regular region of the 
flow. 



Ionian system. We may say that this Hamiltonian system is 
embedded in a larger dynamical system, this time dissipative, 
whose trajectories may or may not converge to zero. If they 
do, the system ends up on the same trajectories as those of 
the smaller embedded system, but in general this need not 
be the case. The general concept of some dynamics embed- 
ded within some other dynamics may be exploited within the 
framework of dynamical-systems theory to design techniques 
that reject some unwanted trajectories of the original dynam- 
ics by making them unstable in the embedding, paralleling 
what the particles do in the fluid-dynamics case. This idea has 
been dubbed bailout embedding of a dynamical system and 
has obvious applications to control theory [Q. 

Trading generality for clarity, in this Section we present 
an example of a bailout embedding for discrete-time dynami- 
cal sytems that closely represents all the qualitative aspects of 
the above described impurity dynamics. Generally speaking, 
given a map Xn+i = T{xn) — x being a point in a space 
of arbitrary dimension — a bailout embedding is the second 
order recurrence 

Xn+2 - T{Xn+i) ^ K(cC„)(cC„+i - T(cC„)), (12) 

where K(a;) is chosen such that |K(a;) I > 1 over the unwanted 
set of orbits, so that they become unstable in the embedding. 
In the discrete system, almost any expression written for the 
ordinary differential equation translates to something close to 
an exponential; in particular, stability eigenvalues have to be 



negative in the ordinary-differential-equation case to represent 
stability, while they have to be smaller than one in absolute 
value in the map case. In order to simulate the dynamics 



of particles, the operator — 
into the particular choice 

K{x) 



(J + 2/3 St"^l) should translate 



-7VT 



(13) 



with 7 = — 2/3St^^,in the map setting. 

To represent qualitatively a chaotic two-dimensional in- 
compressible base flow we choose a classical testbed of 
Hamiltonian systems, the area preserving standard map intro- 
duced by Taylor and Chirikov: 



where 



T : {xn,Vn) (a;„+i,y„+i), 



Xn+i = Xn + — sin(27ry„) 

ZTT 



(14) 



y-n+i 



Vn + Xn+l 



(15) 



and k is the parameter controlling integrability. Recall that in 
general, the dynamics defined by this map present a mixture of 
quasiperiodic motions occurring on the KAM tori and chaotic 
ones, depending on where we choose the initial conditions. As 
the value of k is increased, the region dominated by chaotic 
trajectories pervades more and more of the phase space, ex- 
cept for increasingly small islands of KAM quasiperiodiciy. 
Hence our qualitative description of the impurity dynamics, 
the bailout embedding of the standard map, is given by the 
coupled second-order iterative system defined by Eqs. (p^), 
(HI), (p, and 

Notice that due to the area-preserving property of the stan- 
dard map, the two eigenvalues of the derivative matrix must 
multiply to one. If they are complex, this means that both 
have an absolute value of one, while if they are real, generi- 
cally one of them will be larger than one and the other smaller. 
We can then separate the phase space into elliptic and hyper- 
bolic regions corresponding to each of these two cases. If 
a trajectory of the original map lies entirely on the elliptic 
regions, the overall factor exp(— 7) would damp any small 
perturbation away from it in the embedded system. But for 
chaotic trajectories that inevitably visit some hyperbolic re- 
gions, there exists a value of 7 such that perturbations away 
from a standard-map trajectory are amplified instead of dying 
out in the embedding. As a consequence, trajectories are ex- 
pelled from the chaotic regions finally to settle in the elliptic 
KAM islands. 

To illustrate this. Fig. |] shows the phenomenon in a situ- 
ation in which the nonlinearity parameter has been set to the 
value k = 7. This corresponds to a very chaotic region of 
the standard map, characterized by the existence of minute 
KAM islands within a sea of chaos that covers almost all the 
available phase space. Fig. ||(a) is a close-up — to make it 
visible — of the largest of these islands and the dots there rep- 
resent the successive positions of a set of 1000 fluid parcels — 
evolving according to the standard map — spread initially at 
random over the unit cell. Since none of these parcels were 
initially located inside the island, this is seen as a white spot 
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FIG. 3: The standard map for k = 7 has a chaotic sea covering 
almost the entire torus, except for a tiny period-two KAM island near 
position 0, 2/3. 1000 random initial conditions were chosen, iterated 
for 20000 steps, then the next 1000 iterations are shown. The images 
here are a box -0.05 < a; < 0.05, 0.61 < p < 0.71. (a) Original 
map, (b) 7 = 1.4, (c) 7 = 1.3, (d) 7 = 1.2. 

never visited by the parcels. In contrast, in Figs ^b), ^c), 
and ||(d), the dots are the successive positions, after a number 
of equilibration iterations, of particles initially placed as the 
parcels in Fig. |^(a), but allowing a very small initial discrep- 
ancy 5q = Xi — T{xo) of both dynamics. Notice that now, 
although having started initially outside the island, some of 
the particles settle inside in a process that becomes increas- 
ingly marked as the parameter 7 decreases. 

V. NOISY DYNAMICS 

By virtue of volume preservation, the invariant measure of 
the fluid-parcel dynamics is either uniform, if the system is 
ergodic, or else disintegrates into a foliation of KAM tori and 
ergodic regions, otherwise. In any case, the addition of a small 
amount of white noise, which may be considered to represent 
the effects of small scale turbulence, thermal fluctuations, etc., 
renders the system ergodic with only a uniform invariant mea- 
sure. Thus, the distribution of fluid particles is expected to be 
uniform with or without the presence of noise. The situation is 
however very different if the noise is applied to the dynamics 
of the particles, or, correspondingly, to the bailout embedding. 

Let us consider the foUowing stochastic discrete-time dy- 
namics 

Xn+2-T{xn+i) = VT|,,. (a;„+i - T(a;„)) 

Hn. (16) 



in which, as in the previous Section, x represents the particle 
coordinates and T{x) the fluid parcel evolution. New here is 
the noise term with statistics 

(e«> = 0, 

i^nirn) = e(l - e'^^) 5™.„ I. (17) 

This term forces the particle away from the fluid trajectory 
at every step of the dynamics. However, the actual magni- 
tude of the fluctuations induced in x will be modulated by the 
properties of VT — the flow gradients — along the particle 
trajectory. Notice that for practical reasons we vary the noise 
intensity in correspondence with 7 in order to obtain compa- 
rable fluctuations at different values of this parameter 

Let us again consider the standard map as modelling the ba- 
sic flow, and its noisy bailout embedding that represents a par- 
ticulate impurity subject to both fluid drag and noise forces. 
We are interested in the asymptotic stationary behavior of an 
ensemble of such particles which, invoking the ergodicity of 
the fluctuations, should be well represented by the the his- 
togram of visits that a single particle pays to each bin of the 
space — the full phase space for the basic flow, but a projec- 
tion of the full phase space for the particles — as time goes 

by. 

Figure ^displays a sequence of these histograms in a scaled 
color code for the same nonlinearity parameter fc = 7 as in 
Fig. H, corresponding to the extremely chaotic regime of the 
standard map considered in the previous Section. The se- 
quence of images corresponds to increasing the stability pa- 
rameter exp(— 7). The images make evident the fine filamen- 
tary structure developed by the asymptotically invariant dis- 
tributions due to the combined effect of noise and the ability 
of particles to separate from the basic flow. Remarkably, how- 
ever, these structures appear even in the case in which the 7 
values are larger than those required to produce a spontaneous 
detachment of the particle trajectories. 

Notice that the filamentation here arises from the existence 
of avenues in the phase space that lead to the small KAM is- 
lands on which the particles prefer to stay. A more detailed 
analysis shows that on these avenues the average value of the 
squared separation between particles and fluid trajectories is 
relatively small. Roughly parallel to these avenues, on the 
other hand, there are strips of the phase space that the parti- 
cles avoid. There, the separation between particle and fluid 
trajectories is on average much larger Filamentation is thus 
due to the tendency of the particles to avoid neighboring re- 
gions. 

The same mechanism may also lead to patchiness, not nec- 
essarily filamentary. Looking at flows with weaker chaos, for 
example, for which relatively large KAM islands coexist with 
comparably sized regions of chaos, a situation such as the one 
shown in Fig. |] is typical. There there are relatively small 
avoided regions that separate large patches of larger concen- 
tration of particles around the nonchaotic islands. This picture 
is also testimonial to a property that distinguishes the present 
mechanism from those of Refs. [@, |[|. While there the par- 
ticles group around the unstable manifolds of the homoclinic 
intersections of the basic flows, in our case instead the impu- 
rity dynamics tend to avoid the invariant manifolds. This is 




FIG. 4: Histograms with intensity encoding as the square root of 
invariant probability. Parameters are standard map parameter k — 
7.0, e — 10^^, and, from top left to bottom right, stability parameter 
7 = 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. 



because following these manifolds would mean hitting even- 
tually regions in which the Jacobian eigenvalues are closer to 
one, which locally amplifies the effect of fluctuations on the 
dynamics of the particles. 



VI. DISCUSSION AND CONCLUSIONS 

It has previously been demonstrated that small neutrally 
buoyant particles immersed in a fluid flow, and therefore sub- 
ject to drag forces, may follow trajectories that spontaneously 
separate from those of the fluid parcels in some regions of 
the flow. Specifically, this occurs when the shear is very 
strong compared to the Stokes number associated with the 
particles. Given that in general, the smaller the particles are, 
the stronger is the shear necessary for this phenomenon to 
manifest itself, the actual condition for separation may not be 
fulfilled by some fluid flows of physical interest. However, 
we have shown here that the addition of noise to the forces 



FIG. 5: Histogram with standard map parameter k = 2.0, e — 
10^^, and stability parameter 7 = 0.8. 



acting on the particles may extend the consequences of this 
phenomenon beyond the range of Stokes numbers for which 
separation is possible in an interesting way, namely, the gen- 
eration of inhomogeneities in the asymptotic distributions of 
these particles even in cases where the flow is a highly effi- 
cient mixer. 

There is a large variety of examples in which noise ought 
naturally to be added to the dynamical equations of the parti- 
cles. Thermal fluctuations, for example should be considered 
in a range of small-scale laboratory experiments. The effect 
of small-scale turbulence forcing drifters in oceanographic ap- 
plications might also be another relevant example. In the case 
of plankton dynamics, in which the arisal of inhomogeneous 
distributions is an issue, the autonomous swimming abilities 
of individuals might be viewed as an internal source of noise. 
But could the phenomenon described here be the basis of the 
plankton distribution patchiness? If each individual plankton 
is naively considered as an impurity particle, the answer is ob- 
viously not: their size is far too small for these effects to be 
appreciable. But if there are grounds to consider plankton in 
large-scale colonies moving more or less rigidly in the ocean, 
the consequences of the phenomenon described here need to 
be taken seriously. 

We conclude with an epistemological note: the dynamics 
of neutral particles in flows has suggested to us a generaliza- 
tion to maps that helps to solve some problems in the general 
domain of Hamiltonian dynamics. This generalization in turn 
pays us back by suggesting a way in which inhomogeneous 
distributions may arise in fluid-dynamical problems. We be- 
lieve that this is a remarkable instance OF mutual convenience 
in the interdisciplinary marriage between the two fields. 
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